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A FINITE-STEP METHOD FOR CALCULATION OF THEORETICAL 
LOAD DISTRIBUTIONS FOR ARBITRARY LIFTING- SURFACE 
ARRANGEMENTS AT SUBSONIC SPEEDS 

By James A. Blackwell, Jr. 

Langley Research Center 

SUMMARY 

This report presents a simple theoretical finite- step method for calculating the sub- 
sonic aerodynamic load distributions for arbitrary lifting-surface arrangements. The 
method may be applied to a wide variety of aerodynamic problems involving nonplanar 
wings, wing and pylons or end plates, wing and empennage, and biplanes. The results of 
the theoretical method have been compared with experiment and show good agreement in 
all cases. Applications of the method have been made for a series of illustrative 
examples. 


INTRODUCTION 

Much research is presently being devoted to increasing the range factors of sub- 
sonic cruise airplanes. In order to achieve better range factors through increased aero- 
dynamic efficiency of the airplane design, wings of thinner sections and of increased 
sweep are being used which generally result in additional structural considerations. To 
produce an efficient aerodynamic and structural airplane design, an accurate prediction 
of the aerodynamic forces and their distribution over the wing and other lifting surfaces 
is desirable. 

The results of an accurate theoretical calculation of the lift distribution for a wing 
alone (e.g., ref. 1) are partially invalidated by the presence of a fuselage, pylons, or 
nacelles. The experimental results given in references 2 and 3 indicate the magnitude of 
the interference effect of a vertical surface such as a fence or end plate on the span load 
distribution of a wing. It is shown in these references that the interference effect is sub- 
stantial and should be taken into account. In the airplane design, when the structural 
stresses are required, approximate methods have generally been used to determine the 
lift distributions for interference problems. References 4 and 5 give examples of such 
procedures. Generally, interference solutions are mathematically complex, are restricted 
to particular configurations, and involve much computational labor to obtain an answer of 



satisfactory accuracy. (See, for instance, refs. 6 and 7 in which the calculated interfer- 
ence loads due to end plates and tip tanks, respectively, are shown.) 

It is the purpose of the present report to describe a simple method of lift- 
distribution calculation that is mathematically less complex than existing methods and 
can be applied to arbitrary wings and lifting- surf ace arrangements (i.e., wing and end 
plates, biplanes, wing and pylons, etc.). The method that is used for this analysis follows 
closely the finite- step method used by Campbell (ref. 8) for calculating wing- alone span- 
wise lift distributions. This concept involves using 2N rectangular horseshoe vortices 
placed along the quarter-chord lines of the lifting surfaces and equating the velocity from 
the total vortex system at the three-quarter- chord line to the component of free- stream 
velocity normal to the lifting surface to form 2N equations in 2N unknowns. 

The method is developed in detail and formulas facilitating solution for the span 
load distribution and various other aerodynamic parameters are presented. Applications 
of the method have been made to various lifting- surface combinations with particular 
emphasis placed on a horizontal -wing — vertical-surface configuration. Comparison of 
the results obtained by the present method with experimental data are made. 

SYMBOLS 


A 

b 


aspect ratio, — 
wing span, feet (meters) 


3 H 


■Di 


c L 

c l 


a 


horizontal -tail span, feet (meters) 


total induced drag coefficient, 


Di 




horizontal-surface lift coefficient, 


2S 


horizontal-surface lift-curve slope per radian, 


8Cl 

9a 


C Y 


a 


vertical- surface side-force coefficient, •= ~ 

IpV 2 S 
2 r oo vs 

vertical- surface side-force-curve slope per radian, 


9Cy 

9a 


local streamwise chord, feet (meters) 


2 



average chord, Area of surface feet ( me ters) 
Span of surface 


horizontal- surface section lift coefficient, 

VooC 

dc 7 

horizontal- surface section lift- curve slope per radian, — - 

3 a 

2 r 

vertical- surface section side-force coefficient, 

V TO c 

total induced drag, pounds (newtons) 
efficiency factor 

backwash influence function, foot"* (meter - !) 
defined by equation (6) 

sidewash influence function, foot - ! (meter - !) 
defined by equation (8) 

downwash influence function, foot"! (meter - !) 
defined by equation (11) 

vertical- surface side force, pounds (newtons) 

local deflection at control point in PQ-plane, radians unless otherwise 
indicated 

horizontal-surface lift, pounds (newtons) 
vertical- surface height, feet (meters) 

one-half of number of rectangular horseshoe vortices spread over 
horizontal lifting surface 

free-stream Mach number 
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one-half of total number of rectangular horseshoe vortices spread over 
entire airplane configuration 

coordinate of a particular horseshoe vortex 

coordinate of a particular control point 

coordinates of point on lifting surface in P,Q,R axis system, feet (meters) 

3V& 

area of horizontal lifting surface, feet^ (meters^) 

area of vertical surface, feet^ (meters^) 

semiwidth of horseshoe vortex, feet (meters) 

backwash velocity, in p-direction, feet/second (meters/second) 

free-stream velocity, feet/second (meters/second) 

velocity from vortex system normal to lifting surface at control point, 
feet/second (meters/second) 

sidewash velocity, in q-direction, feet/second (meters/second) 
downwash velocity, in r-direction, feet/second (meters/second) 
longitudinal reference axes 

coordinates of point on lifting surface with respect to given horseshoe vortex 
in X,Y,Z axis system (see eq. (3)), feet (meters) 

OO 


lateral reference axes 



Z,R vertical reference axes 

a local angle of attack at control point in PR-plane, radians unless otherwise 

indicated 


circulation strength, feet^/second (meters^/second) 


r' = 


4ttV 0 


vertical location of lifting- surface root quarter- chord line, measured from 
origin in r-direction, feet (meters) 


V 


spanwise location of vertical surface, measured from origin, feet (meters) 


A 

X 


sweep angle of quarter- chord line, degrees 

taper ratio, Ti P chord , 

Root chord 

longitudinal location of lifting- surface root quarter- chord line, measured 
from origin in p-direction, feet (meters) 

mass density of air, slugs/foot^ (kilograms/meter^) 


0 


angle of lifting surface in QR-plane, degrees 


Subscripts: 


H horizontal tail 

n number designating a particular horseshoe vortex 

vs vertical surface 

W wing 

v number designating a particular control point 
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ANALYSIS 


Basic Concepts 

In order to calculate the subsonic span loading of an arbitrary arrangement of 
lifting surfaces, the lifting surfaces are represented by a system of rectangular horse- 
shoe vortices (fig. 1). For the present analysis, one horseshoe vortex along the chord 
is used. That is, the midpoints of the vortices are placed only at points along the quarter- 
chord lines. An equal number of control points are located along the three -quarter- chord 
lines. 

The velocity from the total vortex system is equated to the component of free- 
stream velocity normal to the lifting surface chord at each control point. Application of 
this tangent-flow boundary condition for a symmetrical loading provides a set of N 
simultaneous equations in the N unknown circulation strengths. Solution of this set of 
equations provides the loading distributions over the lifting surfaces. 

Derivation of Method 

Backwash .- By utilizing the fundamental laws of induced velocity given in refer- 
ence 9 for calculating the induced velocities from a line vortex, the backwash velocity 
induced at a control point Pj, by a singular rectangular horseshoe vortex (fig. 1) on a 
lifting surface of dihedral angle < p may be derived and is given by the expression 

= X F u (x',y,z,s ,0) (1) 


where 


TT ^ _ z cos 0 - y sin 0 f (y + s cos 0) cos 0 + (z + s sin (p ) sin (p 

( x ’)2 + (zcos0-y sin 0) 2 ^ x , } 2 + (y + s cos *)2 + ( z + s sin 0 ) 2 ] V ^ 


(y - s cos (ft) cos 0 + (z - s si n (p ) sin <p 
[(x')2 + (y - s cos 0)2 + (z - s sin 0)2 


and 


x ’ = K - P|T 
y = % - % / 

z = r p~ r n J 


( 2 ) 


( 3 ) 
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The Prandtl-Glauert correction factor is applied to the p- and x- coordinates to 
account for compressibility effects. 

By distributing 2N horseshoe vortices having 2N control points over the lifting 
surfaces (fig. 1), the backwash velocity at any of the control points P^, which results 
from the 2N horseshoe vortices is 

2N 

u (?v’%’ r v) = 4^ 2 Fn Fu ^ n ^ 

n=l 

For symmetrical geometry with respect to the PR-plane, the backwash velocity at 
any of the control points P v becomes 

N 

"(W-) - h 1 r » F u,™ < 5 > 

n=l 

where the function vn is obtained by summing the influence of the vortices on the left 
side of the PR-plane and the influence of the vortices on the right side of the PR-plane, 
that is, 

F u,wi = F u[(p * (%” ^n)’ ( r v ~ r n)> s n> 0nl 

L J Right side 

+ F uT(P^).(V'ln)'( r ^ r n). s n,-d (6) 

L v ' ' ' -i Left side 

Sidewash .- For symmetrical geometry with respect to the PR-plane, the sidewash 
velocity of any of the control points Pj, which results from the 2N horseshoe vortices is 

N 

v (^' w) - h 1 r - K.vn m 

n=l 

with the function F^ vn given by 

Right side 

+ F v|?Pj, - P n ) » (% +( *n) ’ ( r ^“ r n) , s n, “ 0nj (8) 

L -I Left side 



7 



and 


F v (x',y,z,s ,0) = 


x T sin cj) 


(y + s cos 0) cos (p + (z + s sin 0)sin 0 


( x ')2 + (2 cos 0 - y sin 0) 2 < ^,>2 + (y + s cos 0)2 + (z + s sin 0 ) 2] 1/2 


_ (y - s cos <ft)cos <f> + (z - s sin 0)sin 0 ^ 
[<*')* + (y - s cos 0) 2 + (z - s sin 


z - s sin eft j ^ _ x 1 

(y - s cos 0)2 + (z - s sin 0) 2 ( j- (x , )2 + (y . s cos 0) 2 + (z . s sin ^2] 


1/2 


z -t- s sin <fi 


x T 


(y + s cos 0)2 + (z + s sin 0)2 <J Jj x , )2 + (y + s cos 0)2 + (z + s sin 0)2] 1/2 > (9) 


Downwash .- For symmetrical geometry with respect to the PR-plane, the downwash 
velocity at any of the control points P„ which results from the 2N horseshoe vortices is 


w 



r n f: 




( 10 ) 


with the function vn given by 


F w,un = F w[(Pt,-PnHV%)’( r ^ r n)’ s n^n] 


Right side 


+ F w[(Pj,-Pn> {% + %)>( r v- r * )» s n»“ ^ 


Left side 


( 11 ) 
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where 


F w (x',y, z ,s,0) = 


-X 1 COS 0 i yy + » CU5 cp + yz, + a sin tp/sm tp 

(x’) 2 + ( z cos 0 - y sin 0) 2 M( X . ) 2 + (y + s C os 0)2 + (z + s sin 0)2] 1/2 


(y + s cos 0) cos 0 + (z + s sin 0)sin 0 

_ i /< 


(y - s cos 0)cos 0 -j- (z ~* s sin 0)sin 0 J 
[(x’) 2 + (y - s COS 0) 2 + (z - s sin <j>) 2 ] 1/2 


y - s cos 0 a 

(y - s cos 0)2 + (z - s sin 0) 2 < j^, )2 + (y _ g cos 0) 2 + (z . s sin 0) 2] 


172 


y + s cos 0 C V _j 

(y + s cos 0) 2 + (z + s sin 0) 2 < jj x , )2 + (y + s cos 0)2 + (z + s sin 0)2] 1/2 > (12) 


To satisfy the tangent-flow boundary conditions, the velocity from the total vortex 
system normal to the lifting surface is equated to the component of free- stream velocity 
normal to the lifting surface at each control point such that there is no flow through the 
wing at the control points. 

This boundary condition for a wing at dihedral <0 at angle of attack (fig. 2(a)) may 
be expressed as 

Vj_ = -v sin <p + w cos cf> (13) 

For small angles of attack, the tangent-flow boundary condition may be expressed 
as 

v ±(p^> V 1 ^) = V °° sin a v cos $v ~ Voo a v cos <p u (14) 

Defining as and combining equations (7), (10), (13), and (14) yields 

N 

a v = ^ ^n^Wji/n “ ^V.i'n ^ an 
n=l 

Similarly for a vertical-surface deflection i, the velocity from the vortex system 
normal to the vertical surface (fig. 2(b)) may be expressed as 

Vl=-v cos i + u sin i (16) 
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For small deflections and angles of attack, the tangent-flow boundary condition may 
be expressed as 

v ±(P'u’ %’ r v) = V °° cos a v sin V ~ (17) 

Combining equations (4), (7), (16), and (17) yields 

N 

V = £ r n( F i, ra - F v,^n) ( 18 ) 

n=l 

To determine the unknown circulation strengths on a combination of lifting 

surfaces, a set of N simultaneous equations in N unknown circulation strengths as 
given by equations (15) and (18) must be solved. 

For a wing with dihedral <p and a pylon with incidence i, the system of equations 
to be solved is 


M M N 

X ^ = X zL F n( F w,zm " F v,im tan 
v=l v=l n=l 

N N N 

^ = ^ ^ F n( F u,ra ^ “ F v,wi) 

v=M+l i^M+1 n=l 


(19) 


Calculation of Loading Parameters 


The formulas for calculating the airplane loading parameters vary with each air- 
plane configuration. However, as an illustration, the formulas for calculation of the 
loading parameters for a horizontal wing and vertical surface (one on each wing half-span) 
are presented. Out of N control points on the wing half- span and vertical surface, the 
first M of these points are specified to be on the wing half-span and the remainder to be 
on the vertical surface. 


Span loading coeffic ient.- By defining the wing span loading coefficient as 
and using the relation 


r = 


c z cV 


oo 


2 


c Z c /^L c av 

( 20 ) 


the span loading coefficient on the wing may be obtained from the relation 


10 


i ■■■ ■■■ ■ II^MII 


ii-ii-i-ii-mi ■1111111 huh mu ■■■111111 mini hi i i in in i mini m m mi i mu 11 11 i 


I 



( 21 ) 


c 7 c 


kr- 

2 n 


V C L c av/ n 


M 


Z r ^ 2Sn 


n=l 


Similarly, the span loading coefficient on the vertical surfaces may be obtained from the 
relation 




IT, 


n 


\ C Y c avy 


1 N 

n 

> r;2s n 


( 22 ) 


n=M+l 


It should be noted that the circulation or span loading on the vertical surface at the 
vertical-surface — -wing juncture may be determined from the boundary condition that 
specifies the circulation at the vertical- surface root to be equal to the difference between 
the circulation at stations inboard and outboard of the juncture. 

Lif t and side-force coefficients .- The lift of the wing may be expressed as 

pb/2 

L = 2 pVcor(y)dy (23) 

J 0 

Expressing the lift in coefficient form and simplifying gives 

4 pb/2 

c L = v3J 0 r(y)dy <24) 

Introducing r n and changing the integral to a summation results in 

M 

C L = ^ I r n 2s n (25) 

n=l 

Similarly, the side force of the vertical surface may be expressed as 

F y = pV oo r(z)dz ( 26 ) 

^0 


or, in coefficient form, as 


N 

^ 8tt 
C Y = — 

Svs n=M+l 


r n 2s 


n 


(27) 
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DISCUSSION AND APPLICATION OF METHOD 


Using the finite-step method, the spanwise load distribution as calculated for the 
wing alone is compared in figure 3 with the spanwise load distribution of a wing with a 
vertical surface attached to the wing lower surface at the 60-percent wing semispan. The 
calculations show that the addition of the vertical surface to the wing has a substantial 
effect on the wing spanwise loading. 

The magnitude of the interference effects of lifting surfaces is dependent upon the 
geometry of the interfering lifting surfaces. In order to assess the interference effects 
of an arbitrary arrangement of lifting surfaces, equation (19) may be utilized. 

To indicate the capability of the present method for predicting accurately the inter- 
ference effects of arbitrary arrangements of lifting surfaces, comparisons are made with 
experimental data. In addition, several examples are presented to illustrate the range of 
applicability of the method. 


Distribution of Vortices 

The effect of the distribution of vortices spanwise over an arrangement of lifting 
surfaces on the aerodynamic loads as calculated by the present method is shown in fig- 
ure 4 (symbols are used in fig. 4 to indicate distribution). As indicated in figure 4, the 
use of equal span vortices of semiwidth = 0.0250 generally provided a satisfactory 
loading distribution. However, if there are regions of interest for which there are large 
gradients in the span load diagram, it has been found that a closer uniform spacing of 
vortices should be used. In all the following examples, a vortex semiwidth of 0.0250^ is 
used unless otherwise stated. About 1 minute of computing time on the Control Data 6600 
computer system (Langley program A1505) is generally required to obtain a solution when 
50 vortices (2N) are used to describe the configuration. 


Comparison of Theory With Experiment 

A comparison of the calculated results for a 40° sweptback wing of Ayy =4.5 
having a circular end plate with the experimental results as presented in reference 2 is 
shown in figure 5(a). A wing-alone comparison is also made in figure 5(b). Good agree- 
ment is shown with and without end plates. 

The calculated results for a 45° sweptback wing of A w = 2.0 with and without a 
fence are compared with the experimental results of reference 3 in figure 6 for constant 
total lift coefficients. The fence was located at 0.42^ with a height of 0.12^ above and 
below the wing. The agreement is shown to be good. 
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Examples 


Vertical surfaces perpendicular to horizonta l surfaces .- Examples of vertical sur- 
faces perpendicular to horizontal surfaces include pylons, fences, end plates, and vertical 
tails perpendicular to wings. The location, number, sweep, and height of the vertical 
surfaces, in addition to Mach number and vertical- surface cant angle all affect the span- 
wise loading of a given wing. In the following illustrations, unless otherwise stated, the 
chords of the vertical surface are equal to the local wing chord; the wing is flat with 
A w = 30°, A w = 0.33, and A ^ = 6.67; the vertical surface is at zero incidence with 
respect to the flight path and on the wing lower surface. 

The effect of vertical- surface location on the wing spanwise loading is presented in 
figure 7(a). As the vertical surface is moved outboard, the magnitude of the drop in wing 
spanwise loading increases which in turn increases the loading on the vertical surface 
(fig. 7(b)). Since the chords for the vertical surfaces at various spanwise locations vary 
in figures 7 and 8, the parameter is presented instead of Cy to indicate the loading 

on the vertical surface, (it should be noted that 

In figure 8, the effect of the number of vertical surfaces on the wing is shown. The 
drop and magnitude of the wing spanwise loading at the location of the vertical surfaces 
does not appear to be substantially altered by the number of vertical surfaces. Similarly 
the loading on the vertical surfaces remain essentially unchanged. 

The sweep of the vertical surface primarily alters the vertical- surface loading 
(fig. 9). The loading for the sweptforward vertical surface rapidly drops off proceeding 
toward the tip of the vertical surface. This is a result of the sweep removing the tip of 
the vertical surface from the region of the wing influence. 

The effect of vertical- surface height on the wing and vertical- surface loading is 
presented in figure 10. For vertical- surface heights of = 0.0 to 0.2, the magnitude 
of the wing span load distribution appears to be generally proportional to the height of the 
vertical surface. 

By using the Prandtl-Glauert transformation, the effect of Mach number on the wing 
and vertical- surface load distributions is determined and presented in figure 11. 

Parallel liftin g surfaces .- Biplanes, wing and horizontal tail, and wing and canard 
combination are examples of parallel lifting surfaces. As an example of the interference 
effect of parallel lifting surfaces, the effect of a wing on the lift of a horizontal tail is 
shown in figures 12 and 13. The lift coefficient of the horizontal tail is based on the 
horizontal-tail surface area. In figure 12, the lift distribution for a horizontal tail alone 
and a horizontal tail in the wake of a wing is presented. For the surface arrangement 
shown (tail with wing), the horizontal-tail lift coefficient is substantially reduced. The 
magnitude of the reduction in lift coefficient of the horizontal tail due to the interference 
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of the wing agrees with the reduction expected due to the change in tail angle of attack 
resulting from the downwash of the wing (see ref. 10 for empirical formulas for calcu- 
lating the downwash angle at the wing-wake center line). In figure 13, the effect of 
horizontal-tail location behind a wing on the lift distribution of the horizontal tail is pre- 
sented. Moving the tail vertically or further aft of the wing increases the horizontal-tail 
loading. For the horizontal-tail locations shown, there was negligible effect of the hori- 
zontal tail on the wing. 


Optimum surfaces.- It is well known that if the span load distribution of a mono- 
plane wing is elliptical, the induced drag of the monoplane wing will be a minimum, that 
is 


C 


Di 


ireA 


(28) 


with e = 1.00. For lifting- surface arrangements other than the monoplane, it has been 
theoretically shown in references 5, 11, 12, and 13 that efficiency factors greater than 1 
may be obtained and the optimum span load distribution is not necessarily elliptical. For 
example, the optimum wing span loading as given in reference 5 for a 30° sweptback wing 
with a sweptforward vertical surface having 75° of sweep, height = 0.2, and located 

at j^-=0.6 is presented in figure 14(a). The optimum vertical- surface loading for this 

configuration (fig. 14(b)) should be nearly elliptical (ref. 5). Shown also in figure 14 for 
comparison is the span and vertical-surface loading calculated by using the present 
method for the same wing at angle of attack with no twist and the same vertical- surface 
arrangement with zero deflection. Theoretically, e = 1.025 for this configuration with 
optimum loading. In figure 15, the approximate wing incidence and vertical -surface 
deflection as determined from the present method which would be required to achieve the 
optimum load distribution at a Cl = 0.5 are presented. It is interesting to note that for 
this configuration the differential twist outboard of the wing root and just inboard of the 
vertical surface is negligible. This is probably due to the vertical surface acting as an 
end plate, with the result that the area inboard of the vertical surface is more nearly two 
dimensional in character. For sweptforward vertical surfaces, the differential twist 
required to attain an elliptical load distribution on the vertical surface can be seen to be 
substantial (fig. 15(b)); however, for unswept vertical surfaces, probably very little differ- 
ential twist would be required since the induced load is already very nearly elliptical 
(i.e., see fig. 9(b)). 


CONCLUDING REMARKS 

A simple theoretical finite- step method for calculating the subsonic aerodynamic 
load distributions for arbitrary lifting- surface arrangements has been developed. A wide 
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variety of aerodynamic problems involving nonplanar wings, wing and pylons or end 
plates, wing and empennage, and biplanes can be solved. The results of the theoretical 
method have been compared with experiment and show good agreement for the cases con- 
sidered. Applications of the method have been made for a series of illustrative examples 
and good results are obtained. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 22, 1969, 

126-13-01-29-23. 
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(a) Horizontal lifting surface. 






(b) Deflected vertical surface. 

Illustration of tangent-flow boundary condition. 



I 






















0 .1 -2 .3 .4 ,5 .6 .7 .8 .9 


y 

b/2 

(a) Wing* 

Figure 8.- Effect of number of vertical surfaces on wing and vertical-surface loading. = 30°; A w = 6.67; = 0.33; A, 

l = 0.2^: = 0.0; s = 0.0250- 
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